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Analytic Calculation of the Anomalous Exponents in Turbulence: 
Using the Fusion Rules to Flush Out a Small Parameter 
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The main difficulty of statistical theories of fluid turbulence is the lack of an obvious small 
parameter. In this paper we show that the formerly established fusion rules can be employed to 
develop a theory in which Kolmogorov's statistics of 1941 acts as the zero order, or background 
statistics, and the anomalous corrections to the K41 scaling exponents Q n of the nth order structure 
functions can be computed analytically. The crux of the method consists of renormalizing a 4- 
point interaction amplitude on the basis of the fusion rules. The novelty is that this amplitude 
includes a small dimensionless parameter, which is shown to be of the order of the anomaly of 
C2, S2 — (2 — 2/3 ~ 0.03. Higher order interaction amplitudes are shown to be even smaller. The 
corrections to K41 to 0(82) result from standard logarithmically divergent ladder-diagrams in which 
the 4-point interaction acts as a "rung" . The theory allows a calculation of the anomalous exponents 
( n in powers of the small parameter 82- The n-dependence of the scaling exponents £ n stems from 
pure combinatorics of the ladder diagrams. In this paper we calculate the exponents Q n up to 0(<5|). 
Previously derived bridge relations allow a calculation of the anomalous exponents of correlations 
of the dissipation field and of dynamical correlations in terms of the same parameter 82 ■ The actual 
evaluation of the small parameter £2 from first principles requires additional developments that are 
outside the scope of this paper. 



I. INTRODUCTION 
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The aim of this paper is to build on previous work to 
achieve a controlled evaluation of the anomalous expo- 
nents that characterize various correlation and structure 
function in Navier-Stokes turbulence, and in particular 
the exponents Cn that characterize nth order structure 
functions. The main result of this paper is that given a 
single experimental input (for example the value of the 
anomalous exponent of the second order structure func- 
tion), the n dependence of all the other exponents that 
were reliably measured in experiments and simulations 
can be calculated analytically. 

Decades of experimental and theoretical attention (see 
for example |l]-|(|) have been devoted to two types of si- 
multaneous correlation functions; the first type includes 
the structure functions of velocity differences, 



S n (R) = (\u(r + R) - u(r)\ n ) , 



(1) 



where (...) stands for a suitably defined ensemble av- 
erage. A second type of correlations include gradients 
of the velocity field. An important example is the rate 
e(r,f) at which energy is dissipated into heat due to 
viscous damping. This rate is roughly v\Vu(r,t)\ 2 . 
An often-studied simultaneous correlation function of 
e(r, t) = e(r, t) — e is 



K ee (R) = (e(r + R)e(r)) 



(2) 



It has been hypothesized by Kolmogorov in 1941 (K41) 
and 1962 (K62) that statistical objects of this type ex- 
hibit power law dependence on R pLffl: 



In addition, the K41 theory predicted the values of Cn to 
be n/3. Experimental measurements and computer sim- 
ulations show that in some aspects K41 was remarkably 
close to the truth. The major aspect of its predictions, 
that the statistical quantities depend on the length scale 
R as power laws, is corroborated by experiments. On 
the other hand, the predicted exponents seem not to be 
exactly realized. The numerical values of Cn deviate pro- 
gressively from n/3 when n increases §J§. K62 tried 
to improve on this prediction by taking into account the 
fluctuations in the rate of energy dissipation. On the 
basis of a phenomenological model, assuming the distri- 
butions function of energy dissipation to be lognormal, 
K62 reached the predictions 



n nn(n — 3) 
3 18 



(4) 



Besides the fact that these predictions did not follow 
from fluid mechanical considerations, it was pointed out 
H that they are violating basic inequalities that do not 
allow the exponents Cn to decrease, something that al- 
ways happens with (Q) with n large enough. The quest 
for computing the scaling exponents from the equations 
of fluid mechanics was long, arduous, and on the whole 
pretty unsuccessful. 

In this paper we present an approach that is based on 
our own previous findings which culminates in the ana- 
lytic calculation of exponents like C, n and /1. At present 
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the calculation is not completely from first principle. 
We need the input of one number from experiment, say 
82 = C2 — 2/3. Given this number we can calculate all 
the other exponents systematically with 82 being a small 
parameter that organizes our calculations. To first order 
in 82 we recapture ([!]). We will show that the result to 
0(^2) is universal, independent of the details of the cal- 
culations performed below. To second order we find the 
results 

Cn = I - n( "~ 3 W + 28 2 (n - 2)b 2 ] + 0(8 3 2 ) , (5) 
u = 9S 2 (l + 8b 2 S 2 ) + 0(8%) . (6) 

The curves of £„ vs. n are shown in Fig. [I]. We show 
the K41 prediction, the result of our calculation to 1-loop 
order, and the 2-loop result that is presented in Eq. (||). 
While the form of these results is universal, the numer- 
ical value of the dimensionless parameter 62 depends on 
the details of the calculations; wc find that 62 is always 
negative and of the order of unity. Note that the 2-loop 
results correct the unwanted down curving of the 1-loop 
calculation (which is the same malaise as in K62). 
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FIG. 1. The scaling exponents f n as a function of 
n. The calculation is organized by the small parameter 
82 = C2 — 2/3 « 0.03. Shown is the K41 prediction which 
is zero-order in 82, together with our results to first and sec- 
ond order in 62 . To first order the results are the same as the 
phenomenological prediction of K62, and to 2nd order it is 
Eq. (5) with b 2 = -0.55 according to Eq. (72). 

In thinking about the strategy for this work we were 
led by some insights that developed in the context of 
understanding how to compute the scaling exponents of 
the Kraichnan model of passive scalar advection |||| . In 
that model a scalar field T(r, t) is advected by a Gaus- 
sian velocity field u(r, t) which is delta-correlated in time 
but which has a scaling exponent (2 = e. For e = the 
advected scalar has trivial statistics, and for e small the 
model has a natural small parameter. It turned out that 
the calculation of the exponents can proceed along two 
ways. The first, which is non-perturbative, was pioneered 
in pC|]. It considers the differential equations that the 



nth order correlation functions satisfy, and identifies the 
anomalous scaling solutions as the zero modes of these 
differential equations. The calculation of the exponents 
then depends on the calculation of the zero modes them- 
selves, a task that is not at all easy, and therefore such 
calculations were never done for any order but 0(e). In 
this method the renormalization scale is the outer scale L, 
not the inner scale 77, and the dimensionless ratio of scales 
that carries the anomalous part of the exponents of the 
structure functions is L/R where R is defined in Eq. (Q). 
A second method that was discussed in detail in Q con- 
siders instead of the correlation functions the averages of 
higher moments of VT, i.e. (|VT|™). These quantities 
diverge as powers of L/n, and the exponents of this power 
are the same as the anomalous part of the exponent of 
the nth order structure function. The great advantage of 
the second method is that one can write a perturbative 
theory in e for the scaling exponents themselves, without 
any need to compute the zero modes or any other func- 
tions of many variables. Thus the second method allows 
easily computations to 0(e 2 ) Jl2],[l3| and with some more 
efforts to higher orders. The insight gained is that one 
needs to focus on a quantity that offers the most feasible 
calculation, exposing the anomalous part of the expo- 
nents as it appears with a dimensionless ratio of length 
scales. 

In Navier Stokes turbulence the situation is similar. 
On the one hand we have a non-perturbative theory 
which in this case is the infinite hierarchy formed by 
the equations of motion of the correlation function [ p~4[ . 
We can use this hierarchy to demonstrate that anoma- 
lous solutions exist, but the computation of the scaling 
exponents requires a calculation of the correlation func- 
tions themselves. This is a very difficult task that up to 
now was not accomplished in a satisfactory manner. The 
other approach will be described in this paper. It will 
be a perturbative theory for the scaling exponents them- 
selves, not requiring the computation of the correlation 
function along the way. Similarly to the second method 
in the Kraichnan problem it will be based on considering 
limits of correlation functions when p coordinates fuse. 
In that limit we create, even when all the distances are 
in the inertial range, a ratio of large and small length- 
scales that appears raised to the power of the anomalous 
exponent £ n . 

The two previous findings that influence crucially the 
present formulation are the mechanism for anomalous 
scaling that was announced in Ref. Jig] , and the fusion 
rules that were discussed in Ref. |L6]]. In short, Ref. |l5| ] 
exposed the ladder diagrams which appear in the theory 
of turbulence. These diagrams contain logarithmic diver- 
gences that are summable to power laws with anomalous 
exponents. These ladder diagrams contain "rungs" of the 
ladder, that are actually vertices with four, six and more 
"legs", representing 4-point, 6-point and higher order in- 
teraction amplitudes. In Ref. pa] these objects were rep- 
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resented in terms of infinite series of diagrams that could 
not be resummed analytically. This is where the fusion 
rules are now called to save the day. The fusion rules 
determine the asymptotic properties of n-point correla- 
tion functions when subgroups of p coordinates coalesce 
together. As such the fusion rules are non-perturbative, 
and are believed to be exact. We use the fusion rules to 
determine the asymptotic properties of the rungs. This 
is done such that a calculation of (2 from the theory will 
agree with the experimental value of (2- We then show 
that the knowledge of the asymptotics suffices for con- 
structing a calculation of all the other scaling exponents, 
and in particular of C« for n > 2 and of /z. The crux is 
that in the process of determining the analytic form of the 
rungs we discover that their amplitude contains powers 
of a dimensionless small parameter 82 = C2 — 2/3 « 0.03. 
Using this small parameter in the renormalized 4-point 
interaction allows us to develop a systematic expansion 
in orders of 62- At the end of this paper we sketch a way 
to understand the remaining task regarding the origin of 
the small parameter ( 2 ~ 2/3. 

In Sect. U we summarize past results that are neces- 
sary for the present developments. In Sect. Ill we show 



how the fusion rules can be used to determine the proper- 
ties of the rungs in the ladder diagrams appearing in the 
nth order correlation functions. The first important new 
result is demonstrated in Sect. [v| - the numerical coeffi- 
cient contained in the 4-point rung is shown to be small, 
of the order of the anomaly of £2- This result is cru- 
cial since it allows (to our knowledge for the first time) 
to develop a perturbative calculation of the anomalous 
parts of all the other exponents. The physical reason for 



this result is that ( in Sect. tV) we are developing the 



theory around the K41 solution instead of the dissipative 
solution as was always attempted. In Sect. we pave 
the way for the calculation of all the other exponents. 
In Sect. (v| we calculate the scaling exponents by resum- 
ming the logarithmically divergent ladder diagrams up 
to 0(82) (which is known in the field-theoretic jargon as 
the "1-loop order" in the renormalized rungs). We find 
(admittedly to our surprise) that to this order the scaling 
exponents are identical to K62. Like the latter they suffer 
from the violation of the known requirement that Q n can- 
not decrease with n B. In Sect. |\TJ| we show that the 



2-loop order cures the malaise of K62, and we present 
the result ([|) for £„ that in our theoretical estimate is 
valid for n < 12. The exponent // is also computed in 
this Section. If one wanted results for with higher 
values of n one would need to go to 3-loop order (and 
see Sect. VIII where the form of £ n to this order is pre- 
sented), but the present experimental situation does not 
warrant a theoretical prediction of £„ for very high values 
of n. In Sect. VIII we summarize the paper, paying spe- 
cial attention to the range of validity of the theory and to 
demonstrating that no uncontrolled approximations were 
made. 



II. SUMMARY OF PERTINENT PREVIOUS 
RESULTS 

In this Section we present a brief summary of some 
past work which is most pertinent. We refer to Jx3] as 
Paper I, to as Paper II and to jlTj as Paper III. 



A. The basic perturbation theory 

The starting point of the analysis are the Navier-Stokes 
equations for the velocity field of an incompressible fluid 
with kinematic viscosity v which is forced by an external 
force f(r, t): 



-vA' 



u 



P(u-V)u = Pf, 



(7) 



where P is the transverse projection operator P = 
— A _2 V x Vx. It is well known, (and see for exam- 
ple Paper I) that developing a perturbative approach 
]Td| pl[ for the correlation functions and response func- 
tions in terms of the Eulerian velocity u(r,t) results in 
a theory that is plagued with infra-red divergences. On 
the other hand one can transform to new variables, and 
after the transformation (which amounts to infinite par- 
tial resummations in the perturbation theory) one finds 
a renormalized perturbation theory that is finite, with- 
out any divergences in any order of the expansion (cf. 
p2f and Paper I) . One can achieve such a theory using 
Lagrangian variables |23f| ; we find it technically simpler 
to employ the Belinicher-L'vov transformation |22[|, 



v[r \r,t] = u[r + p(r ,t),t)] 



(8) 



where p(r$, t) is the Lagrangian trajectory of a fluid point 
which has started at point r = Tq at time t — to 



p(ro,t) = / u[r + p(r ,T),T]dT 
Jo 



(9) 



The natural variables for a divergence free theory are the 
velocity differences 



w(r \r,t) = v[r \r,t] - v[r \r ,t] 



(10) 



Since the averages of quantities that depend on one time 
only can be computed at t = 0, it follows that the 
average moments of these BL-variables are the struc- 
ture functions of the Eulerian field S n (R) defined by 
Eq. (|l|). It was shown [^2| that these variables satisfy 
the Navier Stokes equations, and that one can develop 
(cf. Paper I) a perturbation theory of the diagrammatic 
type in which the natural quantities are the Green's 
function G a p(ro\r, r', t, t') and the correlation function 
F af3 (r \r,r',t,t'): 
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G af3 (r \r,r',t,t') = 



6{w a (r \r,t)) 



(11) 



f->o 



F af3 (r \r,r',t,t') = (w a (r \r ) t)w p (r \r' ,t')) . (12) 

Physically the Green's function is the mean response of 
the velocity difference to the action of a vanishingly small 
forcing. In stationary turbulence these quantities depend 
on t' — t only, and we can denote this time difference as t. 
The quantities satisfy the well known and exact Dyson 
and Wyld coupled equations. The Dyson equation reads 

d 



[- - vA]G af3 (r \r, r\ t) = G%(r \r, r', 0+)S(t) 



+ J dr 2 G° aS (r \r,r 2 ,0+) J dr x J dh 
x E<5 7 (ro|r 2 ,ri,ti)G 7 ^(ro|ri,r / ,t - h) 



(13) 



where G°^(ro|r, r' , + ) is the bare Green's function de- 
termined by Eq. (3.20) in Paper I. The Wyld equation 
has the form 

oo 

F a f3(r \r,r',t) = J dridr 2 J dtidt 2 G ai (r \r, n, t\) 



D lS (r 1 ~r 2 ,t~t 1 +t 2 ) + $~ f s(r \r 1 ,r 2 ,t-t 1 +t 2 ) 



x G S! 3(r \r',r 2 ,t 2 ) . 



(14) 



In Eq. (13) the "mass operator" S is related to the "eddy 
viscosity" whereas in Eq. ([l4|) the "mass operator" $ is 
the renormalizcd "nonlinear" noise which arises due to 
turbulent excitations. Both these quantities are depen- 
dent on the Green's function and the correlator, and thus 
the equations are coupled. 

The main result of Paper I is a demonstration of the 
property of "locality" in the Dyson and Wyld equations. 
This property means that given a value of \r — ro in 
Eq. (p"3|), the important contribution to the integral on 
the RHS comes from that region where \r\ — t*o| and 
|f*2 — r o| are of the order of \r — ro|. In other words, all 
the integrals converge both in the upper and the lower 
limits. The same is true for the Wyld equation, meaning 
that in the limit of large L and small n these length scales 
disappear from the theory, and there is no natural cutoff 
in the integrals in the perturbative theory. In this case 
one cannot form a dimensionless parameter like L/r or 
r/rj to carry dimensionless corrections to the K41 scaling 
exponents. For r\ <C \r — t*q | <C L scale invariance pre- 
vails, and one finds precisely the K41 scaling exponents: 

G a0 {Xr Q \Xr, Xr', X z t) = X *G a p(r Q \r, r', t), (15) 
F aP {Xr a \Xr, Xr' , X z t) = X^ F al3 (r a \r,r' ,t) . 

One can derive two scaling relations which hold order by 
order, i.e. 



The solution is z — C2 = 2/3. It was also shown that the 
scaling exponent of the Green's function ( [l5| ) is j3 2 = —3. 
Extending such considerations to the higher order struc- 
ture functions leads to the order-by order K41 prediction 
that ( n = n/3. 

Of course, the order by order result (|l5|) which leads to 
( |l6| ) is not necessarily the correct one. If one could resum 
all the diagrammatic expansion one could find nonpertur- 
bative answers that may be different. The whole sum of 
diagrams may diverge when the outer scale goes to infin- 
ity or the inner scale to zero, allowing a renormalization 
scale to creep in even though the order-by-order theory 
is convergent. The difficulty is that no one knows how to 
resum the infinite expansion which exhibits no obvious 
small parameter. 

In this paper we will propose a way out of this dif- 
ficulty. The new thinking is based on the fusion rules. 
Instead of considering fully unfused correlation functions 
only, we will allow some coordinates to be much closer 
together, say within a distance r, whereas the rest will 
be separated by a much larger distance, say of the order 
of R where r <g R. We will show that we can form a 
dimensionless ratio with i?/r, and that such ratios carry 
anomalous exponents that are going to survive the process 
of fusion of coordinates in correlation functions when we 
make structure functions. We will thus be able to recog- 
nize the anomalous exponents even though at first sight 
there is no obvious renormalization scale. 

To clarify how this mechanism works we need to re- 
mind ourselves of the appearance of ladder diagrams in 
the theory of correlations functions. Such diagrams ap- 
pear in the most transparent way in nonlinear Green's 
functions, and we review briefly our past results on these 
objects. 



B. The Nonlinear Green's Functions 

The nonlinear Green's function G 2 ^ 2 (ro\xi,x 2 ,x^,X4) 
describes the response of a product of two velocity fields 
taken at different space-time coordinates to the action of 
two forces /: 



nap-lb, 1 \ / S 2 w a (r \x 1 )w 1 {r Q \x 2 ) 

G 2 2 [r \xi,x 2 , 233,0:4) = ( -rr-, — ; — ztti — i — r 



8fp(ro\x 3 )6fs(r \x4,) 



(17) 



2z + C2 = 2 , z + 2C 2 = 2 . 



(16) 



where for brevity we use the notation Xj = {rj, tj}. Sim- 
ilarly, one defines the nonlinear Green's functions G PjP as 
the response of the product of p velocity differences be- 
tween p distinct points and r to the action of p forces 
in different points. 

In a Gaussian theory (which ours is not) G 2 . 2 would 
be the products of the linear Green's functions like 
G a/3 (r \xi , xs)G' f5 (r \x 2 , X4) . In a non-Gaussian theory 
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it is natural to assume that this quantity is a homoge- 
neous function of its arguments when they are in the 
scaling regime. This means that 



G?J 7 *(ro|Ari, X z h, Ar 2 , X z t 2 , Ar 3 , X z t 3 , Ar 4 , X z t 4 ) 



= X f3i G2 (37S (ro\x 1 ,x 2 ,x 3 ,X4) 



(18) 



From the Gaussian decomposition of this quantity we 
would guess that fa = 2[3 2 = —6. The proof of locality in 
Paper I means that there is no perturbative mechanism to 
change this scaling index. On the other hand, this quan- 
tity, which is a function of four space-time coordinates 
Xi has scaling properties that are not exhausted by the 
overall scaling exponent fa. We have shown in Paper II 
that when we consider its dependence on ratios of space- 
time coordinates in their asymptotic regimes w e pick up 
a set of anomalous scaling exponents. The main result 
of Paper II was that in the regime r± ~ r 2 <C r 3 ~ r 4 
the diagrammatic expansion of this object produces log- 
arithms like In(r3/Vi) to some power. It was explained 
that the sum of such logarithmically large contributions 
is given by (r 3 /ri) A with some anomalous exponent A. 
To make the appearance of anomalous exponents evident 
we review the simplest object that resums to logarithms, 
i.e. the series of "ladder diagrams". 



K *aAA/V * a AA/V 1 

H<h © + 

b AAA/ b AAA/ 
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+ 
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AAA/ *d 



AAA/ d 

AAA/ T©VW * 



+... 




FIG. 2. Diagrams for G 2 ,2(k a , w n , kj,, re c , w c , « d , Wd) 
(Panel (a)) and of the rung R(k a , k;,, « c , (Panel (b)). Di- 
agram (1) in Panel (a) is the Gaussian contribution made of 
a product of two linear Green's functions. Diagram (2) is the 
skeleton contribution, and diagrams (3) and (4) are the 1-loop 
and 2-loop contributions respectively. In Panel (b) we show 
the beginning of the infinite series expansion for the rung, 
with diagram (1) being the bare rung. 




FIG. 3. A diagrammatic representation of the fusion pro- 
cess. A typical n the order correlation function is represented 
by n wavy lines decorating a circle. Panel (a): The fusion of 

2 coordinates to within a distance r which is much smaller 
that the typical separations R between the other coordinates. 
Panel (b): The fusion of three coordinates to within a dis- 
tance r <C R from each other. The availability of all the 
possible diagrams afforded by the topological rules allows us 
to expose a contribution that is precisely G PtP with p = 2 and 

3 respectively. 

The diagrammatic representation of the nonlinear 
Green's function (O) is shown in Fig. ||a, where the no- 
tation of the diagrams is explained shortly in the fig- 
ure legend and at length in Appendix |a|. The 4-point 
"rungs" appearing in the ladders were represented in Pa- 
per II as an infinite series expansion whose beginning is 
shown in Fig. |^b. The rest of the expansion involves ex- 
ceedingly complicated diagrams that we failed to resum 
analytically. Nevertheless, just from general properties 
one could show that the ladder with n rungs contains 
a contribution of order [Aln(r3/ri)] jn\. The summa- 
tion of all these contributions gives a term proportional 
to (r3/ri) A , and this is the observation that we want to 
build on in this paper. 



III. LADDER DIAGRAMS IN A^'TH ORDER 
CORRELATION FUNCTIONS 

In this Section we demonstrate how the anomalous ex- 
ponents of higher order correlation functions can be re- 
lated to resummed ladder diagrams. The idea is to con- 
sider a typical nth order correlation function and to al- 
most fuse p coordinates, p < n, chosen from the available 
n coordinates. The point to observe is that the diagram- 
matic theory allows us to write, upon observation, all the 
topologically possible diagrams appearing in the expan- 
sion of a given object. Thus for example consider Fig. 
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where we represent a general nth order correlation func- 
tion, in which two coordinates are a distance r from each 
other, and all the rest are a distance R from them and 
from each other. While coalescing the two coordinates 
we pull out a contribution that reads Gi,i connected to 
the rest of the diagram with n — 2 coordinates. We know 
that this contribution exists since it is allowed topologi- 
cally. But this fragment is represented in its turn by the 
sum of the ladder diagrams that we present in Fig. @. We 
thus state that whenever we are about to fuse two coordi- 
nates in any nth order correlation function we can expose 
a series of diagrams that are the same as those that ap- 
pear in the expansion of the nonlinear Green's function 
G2, 2(^01^17 X2, %3, £4), together with the logarithmic di- 
vergences that are associated with them. In Appendix |X| 
we explain that in doing so we really take into account 
all the necessary contributions, leaving nothing uncon- 
trolled. 

Now we employ the fusion rules. These tell us that 
as a function of the two fusing coordinates the nth or- 
der correlation function is a homogeneous function whose 
homogeneity exponent is £2- Accordingly, if we succeed 
to resum the ladder diagrams in the limit r <C R we 
should find an anomalous part of £2 from the dimension- 
less power (R/r) A . 

Similarly, we can almost fuse 3, 4 or p coordinates, 
and accordingly pull out of the diagram for the nth or- 
der correlation functions more complex ladders with 3, 
4 or p struts. For example in Fig. ||b we show how the 
fusion of 3 coordinates singles out a fragment that is 
whose ladder diagrams have three struts cf. Fig. |^. The 
fusion rules guarantee that the nth order correlation is 
a homogeneous function of the p fusing coordinates with 
( p being the homogeneity exponent. We will show that 
these more complex ladders also resum into power laws 
in R/r, being responsible for the anomalous parts of £ p . 

At this point all this is a bit formal, since we do not 
have an explicit form of the rungs in the ladder diagrams, 
and we can compute nothing without this knowledge. 
In the next two Sections we will address this issue and 
demonstrate that a judicious use of the fusion rules dic- 
tates enough knowledge of the rungs to take us through 
a useful calculation. 



IV. BUILDING THE THEORY ON THE 
BACKGROUND OF K41 

In this Section we reorganize the theory such that Kol- 
mogorov's 41 theory serves as its "free" limit. In other 
words, we aim at achieving a theory in which resum- 
mations of divergent contributions would directly give 
the anomalous parts of the scaling exponents; the K41 
parts should be obvious order-by-order. This is done in 
two steps, that are correspondingly presented in Sub- 
sects. EVA] and IIVB. 



A. Resummation into K41 propagators 



It was explained in Subsect. II A that our theory is 



developed in the BL-representation, to eliminate spuri- 
ous IR divergences that stem from the sweeping interac- 
tions. The main result of Paper I was that after line- 
resummation each diagram in the BL-diagrammatic ex- 
pansion of the propagators (Green's function and double 
correlation function) converged in the infrared and the ul- 
traviolet regimes. Accordingly, K41 scaling is a solution 
of the order- by-order theory. Nevertheless, the propa- 
gators in the BL-representation lose translational invari- 
ance, and are therefore not diagonal in Fourier space. 
For the purpose of actual calculations it is extremely ad- 
vantageous to rearrange the theory such that the BL- 
propagators become again diagonal in Fourier space. 

The actual resummation that is necessary is presented 
in Appendix [b|. It results in a diagrammatic theory that 
is topologically exactly the same as the standard Wyld 
diagrammatic expansion before line resummation. There 
are two differences as explained in the Appendix [b|. For 
the purposes of our considerations below the main issue 
is the simple form of the the propagators that appear 
as lines in the diagrams: they exhibit K41 scaling expo- 
nents: 



G a/3 (k,uj) = P a p(k)g(k,(j), g(k,uj) = 



F aP (k,u) = P aP (k)f(k,u), f{k,u) 



1 



lu + ij(k) 



,2 1 ~2 



, (is) 



(20) 



In these formulae the scaling exponents are carried by 

7(fc) = c 7 e 1/3 fc 2/3 , <f>(k) = c efc~ 3 , (21) 
where c 7 and are dimensionless constants. 

B. Renormalization to K41 4-point rung 

In this Subsection we determine the form of the 4- 
point rungs of the ladder diagrams in two steps. These 
two steps are based on the following observation: the 
diagrammatic expansion of the rung includes many dia- 
grams, some of which contain in them additional subsets 
of ladder diagrams. In the first step we will consider the 
rungs as if all the diagrams appearing in their infinite se- 
ries were resummed, except for their own internal subsets 
of ladder diagrams. In the second step we will consider 
also the ladder diagrams appearing in the series for the 
rung. We aim to a situation in which all the ladders 
that appear in the theory, like in Fig. ^, contain already 
renormalized rungs. However, instead of evaluating the 
rungs from actual resummations we are going to deter- 
mine their form using the fusion rules. Thus in the first 



G 



step we find the form of the rung that results, upon fu- 
sion, in K41 scaling exponents. In the second step we 
recognize that the rungs themselves have ladders, lead- 
ing to an anomalous correction in the scaling properties 
of the rungs themselves. This being accomplished, we 
will have our final form of the rung. Then we turn to 
the ladder diagrams appearing in the fused correlation 
functions, using the rung as a basic building block of the 
theory. All anomalies of all the measurable statistical 
objects will result from resummations of the remaining 
ladder diagrams. 

Consider Fig. ^|a, in which the rung appears as an ob- 
ject. It is given in terms of an infinite series of diagrams 
in Fig. [2^>. It is in fact a 4-point vertex depending on four 
k- vectors and four frequencies. As a first step we consider 
the value of the rung when all the frequencies are zero, 
denoting it in this limit as R(k a , kb, k c , Kd)- At a later 
point we will explain that this is sufficient for our pur- 
poses. The bare value of this object can be read directly 
from diagram (1) in Fig. |^b, with two bare BL-vertices 
r and one double correlation function. The answer is 



(k a ,kb, k c , Kd) 



S e 



1/3 Y a ^ a (k a , K c , k e )T^^ a (kb, Kd, fe e ) 



(22) 



where fc e = k a — k c = — kb, and Sq is a dimensionless 
constant. 

We demonstrate now that if we use this bare form of 
the rung the fusion rules would predict dissipative expo- 
nents, Cn = n - We first demonstrate this in the context 
of (2. Consider a general n order correlation function 
as in Fig. ||a, and fuse two coordinates, pulling out the 
fragment of 62,2 as shown on the RHS of Fig. ||a. We 
will now compute the scaling exponent by finding the r 
dependence of this fragment when the two coordinates 
approach each other to a small distance r <C R where R 
is the typical distance between all the other coordinates. 
To find the r dependence we must integrate according 
to the explanation in Appendix |^, and to this aim we 
introduce the object T2(r, k): 

Ta(r ' K) = / (£F 4sin2( ^ fca ' r) 

/dto 
-^-G 2 ,2{k ai Ua, -fea, -OJa, «, 0, -K, 0) . (23) 

We are interested in the r dependence of this object in 
the limit re < 1. To calculate T 2 (r,K) in this limit we 
return to Fig. ^|a. Obviously the Gaussian contribution 
diagram (1) is irrelevant in this limit. The skeleton dia- 
gram (2) contributes the following integral: 

xg(k a ,uj a )g(-k a , ~u a )R Q (k a , -k a , k, -k) , (24) 



where the tensor indices of the rung were contracted for 
the longitudinal contribution. The superscript "s" is used 
here and below to denote skeleton contributions. We note 
that in the limit k 3> K the BL vertices are proportional 
to the smallest wavevector k. Thus the rung is propor- 
tional to K 2 /k a 3 ^ 3 - Integrating over the frequencies of the 
two Green's functions g(k a ,u) a ) in this rung [cf. Eqs. (|l§ 
^l| )] results in the evaluation l/[7(fc a )fca 3 ^ 3 ] <x k~ 5 . Thus 
the r dependence of T| is given by 



(25) 



Up to logarithmic corrections this integral is proportional 
to r 2 which is the dissipative solution. Similarly, if we use 
the bare rung in the diagram in Fig. |^t> to determine £3 
we will find £3 = 3. In general we will find = n in- 
stead of the K41 value of n/3. Now one could think that 
the correct values of the inertial range exponents may be 
obtained from resumming all the ladder diagrams with 
the bare rungs. This was the point of view proposed in 
Paper I. In such a case the sought after correction to the 
scaling exponents is of the order of unity, and it is un- 
clear how to develop a controlled resummation. In this 
paper we point out a new way, based on the existence of 
a renormalized rung which gives, upon fusion, K41 expo- 
nents before ladder resummations. The characteristics of 
the renormalized rung in such a scheme are dictated by 
the fusion rules. 

Next we want to determine the renormalized form of 
the rung. To this aim we repeat the exercise of integrat- 
ing over the two Green's functions and the rung, with 
the vertices determined as before in the limit k a m kb 3> 
k c « Kd- But now we leave the exponent of k a in the 
asymptotic evaluation of the rung free, and demand that 

— 5 /3 

the result of the integration will be k a . We find that 
this requires that Rq oc k~ . 

We are now in a position to propose a renormalized 
form of the rung which in a proper calculation could be 
obtained by summing up all the non-ladder diagrams that 
contribute to this rung. This conforms with our basic hy- 
pothesis that all non-ladder diagrams contribute to K41, 
whereas the ladders are responsible for the anomalous 
scaling. Since K41 does not allow L renormalization we 
propose the form 

R ^ (k a , kb, k c , Kd) 

_ £-1/3 -T Q7Cr (fca: K c, k e )T^ Sa (kb, Kd, —fc e ) 



kp k c 



2/3 2/3 
K d 



(26) 



This form gives the K41 overall scaling exponent (which 
is the same as in the bare rung Rq, and in addition 
agrees with the fusion rules for the second order correla- 
tion function with a K41 scaling exponent. In addition 
it is symmetric, as it should be, for exchanging the in- 
dices a and b together with c and d. Note that S is now 
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a renormalized unknown dimensionless parameter which 
will be determined later. 

To proceed, we note that our actual calculation (see 
below) depends really only on the asymptotic proper- 
ties of the rung, which are rigidly determined by the fu- 
sion rules. We can thus attempt to simplify the form 
of the rung as much as possible, preserving the asymp- 
totic and parity properties unchanged. In particular 
we note that the BL vertices T(k a , k c , k e ) have com- 
plicated structure which makes calculations involving 
them rather difficult. Therefore we propose to use in- 
stead Eulerian vertices V(k a , k c , k e ) corrected by a fac- 
tor — 2(fcb • n c )/{k 2 + k\ + ac 2 ]. The correction is aimed 
at reproducing the asymptotic behavior of the BL vertex 
r Q7cr (fc a , k c , k e ) ~ min{fc a , kb, k c }. Thus instead of ( |26[ ) 
one has 

B°&*(k k, k kA - - 4 ^ 1/3 (^-fc e )(K d -fc e ) 

( °' 6 ' c ' d) ~ [*2 + «§ + fee 2 K + «3 + ^] 

V a ^(k a , K C , k e )V^(k b , K d , -k e ) 



fc e 3 («c«d) 2/3 



(27) 



As a further simplification of the actual calculations 
we will use a 1-dimensional reduction of the prob- 
lem (preserving the asymptotic scaling properties and 
parity) in which instead of 3-dimensional integrations 
/ d 3 k/(2iT) 3 we will use 1-dimensional one dk/2-K. 
Then we can disregard the vector indices and replace 
V aicr (k a , k c , k e ) — > k a , (k ■ k') — » kkl (keeping the signs) 
and k 3 in the denominator by \k e \ (because we replaced 
3-d by 1-d integration). The 1-dimcnsional version of the 
rung (p7|) turns into 



R(k a , kb, k c , Kd) 



-4 6 e 1/3 k a k b K c K d \ke 



k 2 ][k 2 



(28) 



Note that here k a , k b , k c , Kd are in the interval ±oo and 
that they carry signs in order to preserve the parity of 
the rungs. 

ka,kb are incoming wave vectors and ac c ,ac^ 
are outgoing, and they conserve momentum, 



ka ~\~ — K c + Kd ■ 



(29) 



Substituting into fl2q) k e — k a — k c — kb — n d one gets 
finally: 



R(k a , kb, AC C , Kd) 

£-1/3 k a kbK c Kd\k a k c \ 



(30) 



'«„ ~ k a K c + Kl][kl - k b K d + K d ]\K c Kd\ 2/3 



In particular, 



R(k, —k, k, — ac) 



-6e^ 3 k 2 K 2 / 3 \k- K \ 
[k 2 — kn + k 2 ] 2 
5 e 1 ' 3 sign(/t ac') k k' \ x / 3 



\k\ 



for ac, k' <C k . 



(31) 



(32) 



To check that we get the right K41 scaling exponents 
with the new renormalized rung we need to recalculate 
the 1-dimensional version of Eq. ( p4[ ) with (^) for the 
rung: 

/°° dk rkr\ 
—4 S m 2 (—)R(k,-k,K,-K) 



— OC 

oo 



duj 
2? 



{k,uj)g(-k, -w) . 



(33) 



Using Eq. (jH 
Thus 



the frequency integral yields —l/2j(k). 



T!(r,K) 



Sk 2/3 roc |fc| 4 /3|AC~Ac|sin 2 



dk- 



[k 2 — kn + ac 2 ] 2 



(34) 



In the limit ac — > the integral simplifies to 



T|(r,K)= 2^/3 /^sin 2 ^), (35) 



where 



5 = 



7TCU 



(36) 



This integral is elementary, reading 

1 , r> /q„ ( 2 N 



T, 



(r, ac) = --J(Acr) 2 / 3 r( - -) « 2~5(Kr) 2 ' 3 , (37) 



where T(x) is the gamma- function. 

The point to notice is that in the asymptotic limit 
ac — ► the only properties of the rung that guaranteed the 
appearance of the scaling exponent 2/3 are the asymp- 
totic properties that we preserved in the series of simpli- 
fications leading to (|30|). In general, we will show below 
that the series of simplifications of the model form of the 
rung are of absolutely no import also for the calculation 
of the anomalous scaling exponents up to 1-loop order. 
We will show below that we get precisely the same ex- 
ponents in this order with any arbitrary analytic form 
of the rung, with tensor indices or without, in 3-d form 
or 1-d form or whatever, as long as the asymptotics are 
preserved, as they are. In 2-loop order this is no longer 
true. The actual numbers obtained in the 2-loop order 
are model dependent. We will show however that the 
sensitivity of the predicted exponents ( n to the model 
for the rung is small as long as n < 8. We need the 2- 
loop order mainly to make sure that it corrects for some 
unacceptable properties of the 1-loop results for higher 
order correlation functions. 
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ka <0) K K 
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AAA/ c 



kc 
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kc 
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+ k b AAA/ 



k c 
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K 
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AAAAA/V- 
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(lc) 
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AAAAAAa- 
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(Id) 



a AAA/ AAAAAAr 



AAA/ 




ic AAAAAA/- 
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FIG. 4. Diagrammatic expansion of G3 : 3 (panel a): Contri- 
butions with no rungs, (diagram 0), one rung and two rungs, 
(panel b): Contributions with three rungs. 



Before we proceed we need to check the self consistency 
of our approach. We need to make sure that all higher 
order non-linear Green's function G p , p (the response of p 
velocities to p forcing) yield, upon fusion, the correct K41 

K41 

exponent for pth order correlation functions £ p = p/3. 
For this aim we have to consider the so called skeleton 
diagrams which are the lowest order connected diagrams 
without loops. For G^^ this is diagram (2) in Fig. ||a. 
for G3,3 the skeleton contribution are shown as diagrams 
(3) on Fig. ||a. Wc must make sure that the skeleton 
diagrams, upon fusion, yield K41 scaling £ p = p/i for 
the appropriate p th -order correlation functions, since our 
grand hypothesis is that the anomalous scaling comes 
only from ladder resummations. 

This test of self consistency is presented in Ap- 
pendix ^|. The important conclusion of this Appendix 
is that the skeleton diagrams for G PyP ({kj, kj}) (with 
asymptotics of the rung defined by the two-point fusion 
rules with £ 2 = 2/3) automatically reproduces the K41 
scaling exponent £ p = p/3 when p points are fused. 



V. SANDING THE FLOOR IN THE ONE-LOOP 
ORDER 

In this Section we demonstrate the most important 
new property of the resummed theory, i.e. that the rungs 
in the ladder diagrams appear with a small parameter. 
This will allow us to develop a controlled ladder resum- 
mation, something that to our knowledge has never been 
available before. In fact, we will lay out in this Section 
all that is needed to calculate the scaling exponents in 
the 1-loop order. In Subsect. VA we demonstrate that 



the prefactor S of the rung is the order of 82 which is the 
anomalous part of £2 and thus small. In Subsect. VB 



we consider the anomalous exponent of the rung itself, 
denoted as 8 ai and stemming from ladder resummations 
withi n th e rung infinite series representation. In Sub- 
sect. VC we reconsider the contribution of the skeleton 



diagrams to the scaling exponents upon fusion, ta king 
into account the anomaly of the rung. In Subsect. VP 



we throw in the inputs: the fact that £3 = 1 and the ex- 
perimental value of (2 ■ The result is Eq. (|6^) which states 
that in the 1-loop order all the unknown parameters are 
numerically identical. From this point the calculation 
of all the other scaling exponents in the 1-loop order is 
straightforward. 



A. The 4-point rung is small! 

Here we show that the coefficient 8 in front of the renor- 
malized 4-point rung ( prj| ) is of the order of the correction 
to K41 of the scaling exponent (2- 



C2 - C 2 



0.03 



(38) 
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To this aim consider the 1-dimensional version of the 
quantity T 2 (r, n) of (p3[) : 



oo 

/dk 1 
^4sin 2 (-fc a r) 



\ / ^-G2,2(k a ,UJa,-k a ,~UJ a ,K,0,~H,0) . (39) 
Z7T 



We will examine the ratio of the contributions of the 
1-loop diagram (3), denoted below as T^(t\k), to the 
contribution of the skeleton diagram Eq. ( |33| ) [diagram 
(2) in Fig. |^a]. After performing the frequency integrals 
the 1-dimensional form of T 2 (r, k) [cf. diagram (3) in 



Fig. Ba] is: 



dk a 2 f 1 , \ 
— tt~t sin -A: a r 

7T7(fc a ) 2 



r 2 (n) (r,.) = l 



8 In 



(— ) 



T 2 (s) (r, K ) 



The sum of all these contributions is as follows: 



r 2 (r, K )=T 2 (s) (r, K )+^r 2 (m) (r, K 



(45) 



(46) 



m— 1 



°° 1 1 

li"(r )K )£if{hi 
m! L VrK 



771 — 



r 2 (B) (r,_«) 

(rft) 5 



We see that, as usual, resummation of the leading con- 
tributions from the logarithmic ladder diagrams results 
in the power function with the exponent 8 which is the 
prefactor of the logarithm in the one-loop diagram, see 
Eq. (fl3|) . Because the expected correction 82 to the K41 

K41 

exponent Q 2 is small ( 82 ~ 0.03 ) we conclude that the 
prefactor of the rung 8 is small as well! This allows us to 
begin with the one-loop approximation in computing the 
higher order scaling exponents £ p with p > 2. 



dq 



27T7(fc a 



-R(k a , -k a , q, -q)R(q, -q, k, -k) . (40) 



In the asymptotic limit defined by nr — > the main con- 
tribution to the integral comes from the 2 symmetric re- 
gions of the q-integration in which \k\ -C \q\ -C \k a \. In 
these regions we use the form (S2I). We calculate 



T 2 (1) (r,«)«25 2 K 2 / 3 



^ Sm ^»)? (41) 



As expected the loop integral over q produces a logarith- 
mic contribution. At this point we use the asymptotical 
identity 



lim 



P°° /k\ /In r°° 

1 / d(kr)f(kr)\n(-) = In ( — ) / d(kr)f(kr 
oJ \kJ \rnJ J 



(42) 



which produces, upon comparison with Eq. ( |35|) the final 
result 



T 2 (1) (r, K ) = 8 In — T 2 w (r, K ) 



1 



(43) 



B. Anomalous correction of the rung asymptotics 

So far we have disregarded the explicit appearance of 
ladder diagrams in the infinite series that defines the rung 
itself. As pointed in Paper II the same kind of ladder 
resummation that is responsible for the anomaly of the 
exponents of the nonlinear Green's functions will also 
contribute an anomalous part to the scaling properties 
of the rung. Nevertheless the outer and inner scale do 
not appear in the rung either, and therefore the anomaly 
is explicit only in the asymptotic regime where we have 
a ratio of large and small scales. In this Subsection we 
flush out this anomaly. 

Instead of (|3^) we expect 

= (47) 
for k 3> ft, K , 

with some anomalous exponent <5 a which is expected (and 
later demonstrated) to be of the order of 8 as it stems 
from the same origin. This correction to the asymptotics 
may be achieved, for example by the following model 
form of the rung ( |3p| ) : 



The factor 8 that was introduced in Eq. (^6|) reappears 
here in front of the logarithm as the effective parameter 
of expansion. 

Analogously one computes the leading contribution of 
the two-loop di agram (4) in Fig. ||a. This is done explic- 
itly in Subsect. VII A and Appendix |d|: 



T 2 (2) (r,ft) = - 



8 In 



(-) 



T 2 (s) (r,ft) 



(44) 



and, in general the leading contribution of the n-loop 
diagram: 



Ra,{ka: kb, ftcj ftci) — R(k a , /Cft, K c K d) | 



|ft C ftd| 



\ka ftc| z 



(48) 



As before, we will argue that the exact analytic form of 
the rung is not important for our calculations, and only 
the asymptotic scaling form is essential. This statement 
will be shown to be exact in the 1-loop order. We thus 
need at this point only to preserve the essential proper- 
ties, i.e. that the outer scale L cannot appear due to 
locality, that the rung has to be symmetric with respect 
a, b — > c, d, etc. 
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C. Contributions of the skeleton diagrams with the 
anomalous 4-point rung 



D. The 2nd and 3rd order correlation functions: 
relations between 5 a , 62 and S 



In this Subsection we reconsider the skeleton diagrams 
appearing in the nonlinear Green's functions G PiP but 
taking into account the anomaly of the rung. In other 
words we are going to compute the scaling exponents 
accounting only for the ladder resummation inside the 
rung, but not the ladder resummation with the anoma- 
lous renormalized rungs. This final step will be done in 
Sects. and [VH. 

We are interested in the scaling exponents of structure 
functions in r representation, and these are obtained from 
the correlation functions in k representation as detailed 
in Appendix Upon fusion we obtain automatically 
contributions behaving like nonliner Green's functions. 
Accordingly the objects of interest in the analysis below 
are the nonlinear Green's functions in which the k de- 
pendence of the fusing coordinates is transformed to r 
representation. The outgoing wavevectors k are left as 
are, and the outgoing frequencis can be put to zero with 
impunity. This results in objects defined in a mixed r, k 
representation, which we denote as T p (r, {Kj}): 

T p (r,{K>})= J Y[^hs(c 1 + ...u Jp )S(k 1 + ...k p ) 



1 (2tt)5 
x /p( r .{ fc i})Gp,p({ fc i> w i.'S.O}) 



(49) 



Here f p {r, {kj}) are one-dimensional versions of the func- 
tions f p (r,{kj}) defined by Eq. jc^). The set {kj} de- 
notes all the outgoing wavevectors. 

We consider the skeleton contributions to the nonlin- 
ear Green's function G PiP , denoted by G s pp . Similarly to 
the definition ( |35| ) we introduce Tp(r, {Kj}) as 

00 p 

T s p {rAk' J })= [ Y[^0^S(u; 1 + ...u Jp )S(k 1 + ...k p ) 

-00 i=1 

X f P (r,{kj})G a PtP ({kj,ujj,Kj,0}) . (50) 

Repeating the calculation of Appendix ^| in the asymp- 
totic regime Kj r <C 1 but with the redefined rung ( fl^ ) 
one gets 



T p {r,{K 3 }) = C p (erYl*r* s *Jl\K j \ l 



/3+5 a 



(51) 



Here C p are dimensionless constants that absorb all the 
numerical factors. In fact, the result ( |5l| ) could be 
guessed directly by recognizing that every rung which is 
connected with the "outgoing" Green's functions G{nj, 0) 



contributes to T p (r, {Kj}) a factor \k. 



l/3+<5 a 



All together 



they give YVj=i l^'l +5a - Convergence of the integrals 
over Kj and ujj implies that neither inner nor outer scales 
may appear, and therefore dimensional consideration re- 
quire a factor r pty 3 +5 *\ 



Consider first the scaling exponent (,2- In Subsect. V A 
we showed that the resummation of the ladder dia- 
grams leads to an anomalous correction to the exponent 
C2 = I which is —5 (cf. Eq. (|46|). But according to 
Eq. ( p 7 !"] ) the ladder in the skeleton contribution brings in 
an additional correction 2(5 a . Altogether we have in the 
1-loop order 



C: 



= §-*• 



2<5 a 



Therefore the exponent 82 defined by (|3 
pressed as follows: 

62 = 2 <L - S . 



(52) 



may be ex- 



(53) 



Another relation between the exponents will follow 
from the analysis of the fusion of three points. To 1- 
loop order the nonlinear Green's function £1*3,3 has the 
skeleton contribution diagram 3 in Fig. ^]a, and the 1- 
loop diagrams in Fig. [|b. The skeleton contribution can 
be read directly from Eq. (pi]): 



(54) 



T$(r,{*j}) = C 3 e>r 3 ^| KlK2K3 r /3+5a 



To discuss the other contributions we refer to Fig. [| in 
which all the diagram of (?3,3 with zero, one and two 
rungs are represented. We have one diagram with no 
rung, three with one, nine with 2. The multiplicity of 3 
in the diagrams of type (2) represent the three possible 
connections of two struts by two rungs. The multiplicity 
of 6 in the diagrams of type (3) represent the different 
pair-permutations of three struts. In general there are 3 n 
diagrams with n rungs, out of which three will have one 
disconnected strut. Diagrams with disconnected struts 
will not contribute in the asymptotic regime that inter- 
ests us here. Thus out of the diagrams in Fig. ||a only the 
skeleton diagram (3) remains in the asymptotic regime. 

In general, with n rungs we have 3™ — 3 fully linked 
diagrams. This number is 6[(3™~ 1 - l)/2], and the num- 
ber [(3" _1 — l)/2] counts the topologically distinct fully 
linked diagrams with n rungs. Thus for example we rep- 
resent in Fig. ^> the four topologically distinct contri- 
butions with three rungs. These are all the 1-loop lad- 
der diagrams contributing to the 3'rd order correlation 
function. We show now that of these four terms dia- 
gram (la) does not contribute a logarithmic divergence, 
whereas the other three contribute the same logarithmic 
term. In fact, this is the beginning of a systematic rule: 
the only diagrams that contribute logarithmic terms in 
the 1-loop order are those in which the last rung appears 
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to the right of the skeleton diagrams. Similar rules will 
be established below for higher loop contributions. 

Consider then the 1-loop diagram (la) in Fig. |]b in 
which this rule is not obeyed. We focus on the loop made 
by the two rungs and the Green's functions q and q' ', 
considering the asymptotic regime k a ,kb >• q ^> K a ,Kb- 
This is the only regime in which a logarithmic divergence 
is possible. In this regime q' ~ q" w k a + kf,. Thus 
the rung R a (k a , fc&, q, q') contributes to the loop g 1 / 3+l5 ». 
The rung Rq,q', Ka ,q" and the Green's function G(q') do 
not contribute any q dependence to the integrand. The 
uj integration over the product of the Green's functions 
G(q) and G(q') gives approximately l/j(q') and again 
contributes no q dependence. Finally we have the evalu- 
ation 



T 3 (la) oc 



dqq 



l/3+Sa 



(55) 



Clearly, this diagram does not exhibit a logarithmic di- 
vergence and as such it does not contribute to the renor- 
malization of the scaling exponent. 

The other three diagrams in Fig. |]b (namely lb,lc and 
Id) are different, they all have a logarithmic divergence. 
The reason for the difference is that in these three dia- 
grams there are four Green's functions, instead of five in 
diagram (la), which carry large wavevectors. This is the 
same situation as in the skeleton diagram (3) in Fig. |^a. 
In the loop we have now two Green's functions, instead 
of one in diagram (la), that carry small wavevectors q. 
This difference leads to a different q dependence in the 
loop, and to a logarithmic divergence. We demonstrate 
this explicitly in the next paragraph, but we already draw 
the conclusion which is general: 1-loop ladder diagrams 
with logarithmic divergences are those in which the ad- 
ditional rung (compared to skeleton diagram) has been 
positioned to the right of the skeleton structure. 

Explicitly, consider diagram (lb) in Fig. ||b. The 
rung R a (k a ,kb,q,q') contributes g 1 / 34 ^ as before. But 
now also the rung R a (q' , k c , —q, n c ) contributes the 
same g-dependence. On the other hand the rung 
R a {q, —q, K a , Kb) contributes |g|~ 1 ~ 2l5a . The u> inte- 
gration with the product of the two Green's func- 
tions G(q, w)G(—q, —us) is the same as (|33| ) leading to 
—1/27(3). 111 total we have a logarithmic integral. 

The diagram (Id) is very similar to (lb); it has the 
same rung structure at the left, and the rightmost rung 
is R a {q, —q, i^ai K c)- This makes no difference to the q de- 
pendence and thus to the logarithmic divergence or to the 
factor in front of the logarithm. Diagram (lc) is slightly 
different, having the third rung on the same ladder as 
the second rung. Nevertheless the rung R a {q' ', k c , q, —q)) 
contributes exactly the same q dependence as the two 
rungs in diagrams (lb) or (Id). Thus it yields at the 
end the same factor with the same logarithm. Finally, 
comparing lc to diagram 3 in Fig. ^|a we see that the 
loop structures are identical in both, and thus if diagram 



3 had a prefactor 5, we can immediately conclude that 
the three diagrams (lb, lc and Id) will result in a total 
prefactor of 35: 

1 



TP(r, { Kj })= 36 In — T^(r, {«,•}), 



(56) 



where k = [K1M2K3] • The leading contribution from 
the higher loop diagrams can be seen to contribute higher 
order terms in the series of a power law, similarly to the 
mechanism displayed in Eqs. (p^j-|4rj|) : 



T 3{r,{Kj}) = 



Tt\r,{^}) 



Substituting ( |54| ) we find finally 

r 3 (r,{ Kj }) = C 3 e-H 1+3 ^-^ 
Accordingly to 1-loop order we write 
C 3 = 1 + 3S a - 36 . 



(57) 



(58) 



(59) 



At this point we use the exact, nonperturbative result 
that £3 = 1 to find the relationship between S a and 
5: S a = 5. Together with (53) we get the important 
conclusion that all our 8 , s are the same: 



62 = 5 a = 5 = £2 - C; 



K41 

2 



0.03 



(60) 



We should stress that this important result is obtained 
using only the asymptotic scaling properties of the rung. 
Changing the explicit form of the rung without ruining 
the asymptotics will affect only the subleading terms in 
the analysis. The leading logarithmic terms are insensi- 
tive to the details of the analytic form of the rung. 



'AAA/ T>WWW 

2 AAA/ UvW 

3 AAAAAAAA7V 




AAA/ 



AAA/ 



4 AAAAA/V- 



AAAAAA/- 



1 AAAAA7V 



2 AAA/ AAA' 



+ 




3 AAA UAAAAAAAAAr- 



4 AAAAAAAA/V 



+ all permutations of ( 1,2,3,4) lines 



FIG. 5. The skeleton contributions in the diagrammatic 
expansion of G4.4. 
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VI. ANOMALOUS SCALING EXPONENTS IN 
THE 1-LOOP APPROXIMATION: SURPRISE, 
SURPRISE 

We are poised to compute now the anomalous correc- 
tions to all the scaling exponents of the p-order correla- 
tion functions in the 1-loop approximation. Start with 
the 4th order nonlinear Green's function, and consider 
the skeleton diagrams in Fig. ||. In the one loop order, to 
obtain a logarithmic divergence in the asymptotic regime 
we must add the additional rung on the right of the skele- 
ton structure. The combinatorics are elementary: Each 
skeleton diagram can host a new rung on the right in six 
different ways. Once a rung has been put in place the 
leading (logarithmic) contribution to the loop integral is 
the same as the loop integrals considered in the last Sec- 
tion. It gives the same logarithm with the same prefactor. 
The only difference is in the combinatorics. We can thus 
write by inspection 



T^(r,{ Kj })^65\n — T 4 w (r, {«,•}) , 



(61) 



where n is the geometric mean of all the Kj. Resum- 
ming the leading contributions of the higher order loop 
diagrams results in the power law 



^(r, {/% }) 



T 4 (s) (r,{^}) 



16(5 



(62) 



Using Eq. (|5l| ) we reach the final result 

d = 4/3 + AS a -65 = 4/3 - 25 2 , 1-loop order. (63) 

The analysis of the 1-loop order contribution to the 
anomalous exponents of the p-order correlation functions 
is as straightforward. There are p(p — l)/2 possibilities 
to append an additional rung to the right of the skeleton 
structure of the p-order nonlinear Green's function. All 
these diagrams contribute identical leading order loga- 
rithmic terms, with the same prefactor, summing up to 
an anomalous correction to the scaling exponent of the 
skeleton diagrams which is 5p(p — l)/2. According to 
Eq. (|5l|) the scaling exponent of the skeleton contribu- 
tion itself is corrected with respect to K41 by p5 a - Thus 
altogether 



P 



(p = ~ + P$a - 5 



(64) 



Using Eq. <M) 



P P\P ) 

( p = - - 5 2 , 1-loop order. (65) 

We note that this formula, which is valid in our case 
to 1-loop order only, is identical in prediction to Kol- 
mogorov's log-normal phenomenological model (known 



as K62). This is interesting, as it stems from the non- 
trivial topology of the ladder diagrams, in which only 
the most leading were considered. The present authors 
find the connection between lognormality and ladder di- 
agrams unexpected. 

Nevertheless we should recognize that in the present 
approach this result has a limited region of validity. The 
analysis of the 2-loop order which is provided below will 
show that Eq. ([|) is only valid when p5 2 <C 1. The 2- 
loop order will contribute positive terms of the order of 
5\p 2 (j> — 3), reducing the negative tendency of the correc- 
tion to K41. Accordingly the present theory will not suf- 
fer from the well know deficiencies of the K62 log-normal 
model which for us is only a first order result. 



VII. ANOMALOUS SCALING EXPONENTS IN 
THE 2-LOOP APPROXIMATION: K62 IS CURED 



In this Section we calculate the 2-loop contributions 
to the scaling exponents £ p . Even though these contribu- 
tions are very small when p5 2 is small (for, say, p < 6), 
they become important for larger values of p where K62 
begins to turn down the n dependence of £„ . In addition 
this calculation allows to present clear ranges of validity 
for the 1-loop and 2-loop calculations. 



A. 2-loop contributions to £2 

We consider the 2-loop diagram (4) in Fig. ||a. Sub- 
stituting it instead of G 2j2 in Eq. ( |39| ) we obtain the 

quantity T 2 ( r J K )- We want to compute the correction 
that this diagram gives to the skeleton diagram (2), and 
to this aim we divide it by TJ (r, k) . In the asymptotic 
regime nr <C 1 the loop integrals over q\ and q 2 con- 
tribute mostly in the range k ^> qi,q 2 S> k. In this 

regime the integrals over k a , w a cancel from the ratio 

( 21 

of T 2 V J (r,K)/T^(r,K). In addition the Green's functions 
G{k c ) and G(nd) cancel. Thus this ratio can be read 
from the ratio of the corresponding diagrams for 6*2,2, 
or taking the diagram (4) and amputating the incoming 
and outgoing legs. We still need to divide by the rung in 
diagram (2) where k a is replaced by 1/r: 

T 2 (2) (r, k) _ 1 f dq 1 dq 2 J(qi) J{q 2 ) 

7?(r,/s) ~ R(l/r, -1/r, k, -k) J {2Tr) qi q 2 

— OO 

Xij(r _1 ,-r _1 , qi-qi)R(qi,-qi,q 2 -q 2 )R(q 2 -q 2 , K-k), 

/CO J I 
^L G {q,u)G{- q -u>) = -^ ) . (66) 

In Appendix ^ we analyze this integral in the asymptotic 
limit sr <C 1 with the final result 
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{-) 


+ 61 




\Kr/ 





where b\ is a dimensionless constant 
61 « -0.434 . 



I5(r,«) 



(67) 



(68) 



In Eq. (^) the In 2 term accounts for the exponentiation 
of the 1-loop contribution, whereas the In term provides 
the 2-loop correction to the scaling exponent ( 2 ■ Instead 
of Eq. (p3J) we now read 



6 2 = 2S a - 6 - hS 2 



(69) 



A second relation between these exponents will be de- 
rived in the next subsection. 



B. 2-loop contributions to £3 

The calculation to 0(5 2 ) of the contributions to £3 and 
of higher order £ p due to ladder resummations introduces 
for the first time 6-point irreducible interactions ampli- 
tudes. These appear in the ladder diagrams as rungs with 
six legs, arising from diagrams that due to their topology 
cannot be resummed into reducible contributions consist- 
ing of two 4-point rungs and one Green's function. The 
6-point rung is discussed in Appendices [A| and ^| In par- 
ticular in Appendix |G] we explain why the functional de- 
pendence of ( p on p can be understood completely on the 
basis of the analysis of ladders with 4-point rungs. This 
stems from the fact that the reducible and irreducible 
contributions to the 6-point rung are of the same order, 
and their combinatorical factors are identical. 

There are many possible two loop diagrams involving 
4-point rungs that appear in the expansion of G^^. How- 
ever, we are only interested in those contributing a loga- 
rithmic divergence in the asymptotic regime. As before, 
to get the relevant diagrams we need to append the last 
rung to the right of the existing 1-loop structure. Thus, 
we begin with the 3 logarithmic diagrams in Fig. ||b (i.e. 
lb, lc and Id) and consider all the diagrams that are ob- 
tained by adding an additional rung on the right which 
connects two struts. The nine resulting diagrams are 
shown in Appendix These diagrams are subdivided 
into two groups: three diagrams in which the last rung 
connects the same struts as the previous rung, and six di- 
agrams in which the last rung connects different struts. 
In Appendix [f] we explain that the first group of dia- 
grams gives exactly the same asymptotic integral as the 
2-loop contribution to £2- This statement should be reit- 
erated because of its importance to the structure of the 
theory: the integrals are different, but once the limits 
k ^> qi , q 2 3> k are taken, the resulting integrals coalesce 
with those computed in the previous Subsection. Thus 
the contribution to C3 from these three diagrams will be 



3b\S 2 (which is 3 times larger than the corresponding 
contribution to £ 2 ). 

The six diagrams of the second group look topologi- 
cal^ different, but again in the asymptotic regime coa- 



lesce into an identical integral Eq. (El) with VP — > v]/ 2 . 
where 



*2 (01,(62) 



(70) 



<Ziki| 1/3 |<7i 



g 2 | 4/3 sign(g 2 ) 



(|<?l| 2 / 3 + |32| 2 / 3 + I<?1 + q2\ 2 / 3 M + 9291 + q 2 ) ' 

Following the procedure outlined in Appendix |e] we find 
the coefficients of expansion 

a 2 = 1 , b 2 w -0.55 . (71) 

Finally we get the 2-loop form of £3: 

Ca = l + 3S a -3S~5 2 (3b 1 + 6b 2 ) . (72) 

Demanding again £3 = 1 we find from Eqs. (^, Jn|): 



6 = S 2 - (h + 4b 2 )6 2 + 0(6 3 2 ), 
S a = S[l + (61 + 2b 2 )5] . 



(73) 
(74) 



These results are used in the next Subsection to calculate 
( n for n > 3 to 2-loop order. 



C. 2-loop contributions to £ p , p > 4 

The calculation of the contribution of 4-point rungs to 
C p for higher values of p does not necessitate the evalu- 
ation of new integrals. In Appendix |f] we explain that 
all the 2-loop integrals appearing in the ladders of Ga^ 
and higher order nonlinear Green's functions are identi- 
cal in the asymptotic regime to one of the two integrals 
appearing in the 3-order quantity. The only differences 
are in the combinatorial factors that account for how 
many ways we can choose the rungs to connect between 
p struts. 

If the second rung is connecting the same struts as the 
rung before it we have the same combinatorial factor as 
in the 1-loop order, namely p(p — l)/2. This provides a 
contribution to ( p which is p(p — l)<5 2 6i/2. If the second 
rung is not connecting the same struts as the rung before 
it, we have p(p — l)(p — 2) contributions. This is due to 
the existence of p(p — l)/2 ways connect two struts with 
the first rung, and then 2(p — 2) ways to connect one of 
these two struts with the remaining (p — 2) struts. This 
leads to a contribution p(p — l)(p — 2)b 2 S 2 . We should 
stress that the loops must have a joint strut to give a In 
contribution. Two disconnected loops lead only to In 2 
contributions, which do not enter the 2-loop corrections 
to the scaling exponents. In total we find 
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-p(p-l)( P -2)b 2 S 2 + 0(S 3 ) . (75) 
Substituting Eqs. ( f73| , f74|) we obtain finally 

<p = f " ^^[l + 25 2 6 2 (p - 2)] + 0(<5 2 3 ) (76) 

We should stress that the functional form presented in 
this equation is solid. It is shown in Appendix ^| that 
the contribution coming from 6-point irreducible rungs 
is only renormalizing the value of b 2 which anyway de- 
pends on the precise analytic form of the 4-point rung 
which is not available at the present time. We estimate 
the range of validity of this order of the calculation by the 
Holder inequalities, which disallow a nonlinear increase 
in the Cp as a function of p. The inflection point where 
this requirement is violated may serve as a good estimate 
for the range of validity. This inflection point occurs at 
p sa 1.4 — 1/(6 62 #2) ~ 12. In Fig. |l| we show, within 
this range, the K41 prediction, the 1-loop approximation 
(equivalent to K62) and our 2-loop final result. It is ob- 
vious that the 2-loop loop prediction goes considerably 
beyond the range of validity of the K62 formula which has 
an unphysical maximum at p ~ 11.5. We believe that all 
the reliably measured values of £ p agree very well with 
this prediction. 

Using the bridge relation H [i = 2 — ( 6 we predict 



f i = 96 2 (l + 8b 2 6 2 ) 



(77) 



Plugging in the numbers we get fi — 0.235 + 0(<5 2 ). This 
is to be contrasted with the K62 prediction fi w 0.27. We 
conclude that the 2-loop contribution is very significant 
for experimentally measured exponents. If one wishes 
to obtain theoretical results for £ p with higher values of 
p one needs to consider the 3-loop contributions, which 
pose no further conceptual difficulties. Nevertheless the 
experimental situation does not warrant at the present 
time the effort needed to accomplish such a calculation. 

One should stress before closing this Section that the 
form of Eq. (76) is universal, stemming from the struc- 
ture of the ladder diagrams and from combinatorics only. 
However the numerical value of b 2 is model dependent. 
We have checked that changing the form of the rung keep- 
ing the asymptotics unchanged results in b 2 remaining 
negative while its value not changing by more than a fac- 
tor of 2 or so. At this moment in time one can determine 
b 2 using the value of (4 from experiments, allowing us 
then to predict accurate values of £„ for n up to 12. It is 
our plan however to develop in the near future a theoret- 
ical equations for the 4-point and 6-point rungs, leading 
to an ab-inito determination of their analytic forms, and 
with them of the parameters S and b 2 . 



VIII. SUMMARY AND DISCUSSION 

The main steps of this and previous papers leading to 
the present results have been as follows: 

• The theory is developed using BL- velocities to elim- 
inate the spurious infrared divergences that are due 
to sweeping effects when Eulerian velocities are em- 
ployed. 

• The Dyson- Wyld perturbation theory was line re- 
summed in order to achieve order by order conver- 
gent perturbation theory with K41 propagators as 
the lines in the theory. At this point the objects 
of the theory are two 2-point propagators (Green's 
function and correlator) and one 3-point vertex. 
The 3-point vertex is in no way "small" , and renor- 
malizing it does not change this fact 

• Multipoint correlation functions are considered 
when p coordinates coalesce together. In the fu- 
sion limit kt — > it is advantageous to reorganize 
the theory in terms of one propagator (K41 Green's 
function), and 4-point, 6-point vertices etc. (the 
rungs) . The series of diagrams contributing to the 
fusion limit are then simple ladder diagrams. 

• The crucial step of the theory is achieved by two 
requirements: (i) the 4-point rung should be consis- 
tent at the level of the skeleton diagrams with the 
fusion rules with K41 scaling exponents, (ii) The 
resummation of the ladder diagrams that appear 
when 2 coordinates fuse together should lead to the 
correct value of £ 2 . These double requirements ac- 
complish two things in one go: (i) the theory is now 
developed around the K^l limit, leading to the ap- 
pearance of the small parameter 8 2 in front of the 
4-point rung, and (ii) all the anomalies are coming 
from the ladder resummations. The 6-point rung 
is shown explicitly (Appendix |g|) to be of second 
order in the small parameter, 8-point rungs are of 
third order, etc. 

• We computed the anomalous exponents in 1-loop 
order, inputting the value of £ 2 and requiring that 
£3 = 1. The result is that the scaling exponents are 
predicted to this order to agree with the log normal 
model K62. We showed that to this order the result 
is universal, independent of the simplifications and 
of the model form of the rung. 

• We computed the anomalous exponents in 2-loop 
order. The malaise of K62 is cured, the 2-loop con- 
tribution has a sign that lifts up the exponents from 
the down curve of the K62 parabola. While the 
form of the 2-loop result is universal, the numerical 
value of the parameter b 2 appearing in the final re- 
sult is model dependent, with contributions for the 
4-point and 6-point rungs. 
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To improve upon the present theory one needs to de- 
velop a theory for the 4-point and 6-point interaction 
amplitudes. Here we determined only the asymptotic 
properties of the 4-point rung, and this allowed us to 
predict the form of the scaling exponents, but an input 
of the value of the anomalous part of £2 was needed to 
nail the 1-loop order. In fact we could use the value of 
£4 to fix the value of b 2 and gain a solid prediction of all 
the exponents to 2-loop order. Such a prediction for Q n 
would be valid up to n « 12. It is very easy to generalize 
the result that we have to 3-loop order, with the intro- 
duction of yet one more parameter associated with the 
3-loop integrals, say 63, which included also contributions 
from the irreducible 8-point rung. The result would read 



n n{n - 3) 

Cn= 3 T~ h[1 



6S2 &3 (n 



l)(n 



25 2 (n-2)b 2 
2)] + 0(51). 



(78) 



We stress that this form stems from the structure of the 
ladder diagrams, and we consider it very solid. From one 
point of view we can now use the value of £5 to fix 63 to 
provide a prediction that is valid for any n within exper- 
imental reach for quite some time. But this is not the 
main point. The main point is that we have identified 
the coefficients appearing in this formula with particular 
objects, i.e the 4-point and higher order vertices which 
appear in the theory as the rungs of the ladders. Obvi- 
ously, a calculation of the renormalized rungs from first 
principle would remove the need to input experimental 
information altogether, affording us a complete theory 
of the scaling exponents of isotropic turbulence. At this 
point this is still not in the cards. 
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APPENDIX A: EXPLANATION OF THE 
DIAGRAMMATIC EXPANSION IN FIGS. 2, 3 

It is important to stress that in the present theory we 
take into account all the necessary contributions. To un- 
derstand this we need to say a few more words about the 
representation of G2,2 in terms of simple ladders only, 
and of the nth order correlation function with p fused 
coordinates as shown in Fig. |^. 
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FIG. 6. The contributions to the fusion of 2 and 3 coordi- 
nates from diagrams containing two and three correlators in 
the cross section. 

The presence of these diagrams also complicates the 
discussion of the fusion rules and Fig. [|. When we pool 
out the fragment containing the two coalescing coordi- 
nates we will generate in the series expansion of this frag- 
ment also the diagrams containing two correlators in the 
cross section. To understand their roles we will always 
consider the fragment of the diagrams to the right of this 
special cross section as belonging to the main body, see 
Fig. H But now the fragment pulled out becomes a 4th 
order correlation function instead of (j2,2- One can ana- 
lyze the role of this diagram in the fusion limit, and the 
conclusion is that when r <C R it contributes the same 
power (R/r)^ as the diagrams considered in the body 
of the paper. We thus conclude that the procedure fol- 
lowed in the body of the paper is amply sufficient for the 
calculation of Q 2 ■ 

The same type of considerations apply when we fuse 
p coordinates. In the series expansion for G PiP we will 
have however ladder diagrams whose topology requires 
resummation into 2p-point irreducible interaction ampli- 
tudes which serve as new types of rungs. The first one 
appears at the level of G^, and is the 6-point rung that 
is discussed explicitly in Appendix ^. It will be argued 
that the 6-point rung is of second order in the smallncss 
S 2 that characterizes the 4-point rung. Respectively the 
2p-point rungs are of 0{5 2 ~ 1 ) and thus rungs with p > 3 
will not affect the analysis of this paper. 

On top of the diagrams considered in the body of the 
paper one needs to consider those having p correlators in 
the cross section. The conclusion is however the same - 
the diagrams considered give all the necessary informa- 
tion, no new (or more divergent) information is available 
in the diagrams that we do not consider explicitly. 

At the end of the day the procedure of calculating 
the scaling exponents to 0(5 2 ) leaves us with only three 
dressed objects in the theory, the dressed 2-point Greens' 
function and new 4-point and 6-point dressed vertices 
which we call the rungs. We will show that the rungs 
are small and therefore these three objects suffice for a 
consistent and controlled theory of the scaling exponents. 

A reader who is an expert in diagrammatic theo- 
ries may feel worried that the procedure advocated here 
mixes into the 4-point vertex some 2-particle reducible di- 
agrams. This is indeed so. In other words (and see Paper 
III for more details) there are contributions in which the 
cross section includes two propagators, like one Green's 
function and one correlator, or two Green's functions ori- 
ented oppositely to the ones plotted in Figs. || and [|, 
which we resum into our 4-point vertex. Usually one 
would prefer not to include 2-particle reducible diagrams 
in the 4-point vertex, but rather to distinguish differ- 
ent 4-point vertices. In our case one would need three 



different 4-point rungs, all presented in detail in Paper 
III. Indeed, if we planned to compute our 4-point vertices 
from summing up diagrams it would be very advisable to 
distinguish different types of rungs. In this paper we de- 
termine however the properties of the 4-point vertex from 
the fusion rules, and it makes no difference at all how we 
classify the diagrams. For the sake of clear presentation 
it is much better to have just one type of rung. This rung 
is 2-particle reducible and hides three types of 2-particle 
irreducible rungs as a single object. The conclusions of 
the analysis are independent of this simplification. 



APPENDIX B: RESUMMATION INTO 
DIAGONAL K41 PROPAGATORS 

The starting point of this rearrangement are the mass 
operators in k,u> representation Ti a p(ro\k-\_, k 2 , w) and 
^apirolki, k 2 , u>). Define the "diagonal" part of the mass 
operators as 



fci + k 2 

CaM ^ ) 



<?W( 



2 ' J (2tt) 
fci + fc 2 . _ f d(ki - k 2 ) 



d{kl 3 fc2) S a/3 (fci,fc 2 ,0), (Bl) 



-) 



(2tt) 2 



*a/3(fel,fc 2 ,0) . (B2) 



In these definitions Tq disappears. The reason is that 
for objects which are time independent the Eulerian and 
BL-representations are equivalent and the designation Tq 
is unneeded. Here we have objects with u) = 0, or time- 
integrated quantities. It was shown in ref. (2f| that time 
integrated quantities are related to simultaneous correla- 
tions, and as such they lose the r designation. 
Denote the rest of the mass operators as 

S Q/ 3(ro|fci,fc 2 ,^) = Sa/?(ro|fei, k 2 ,u>)a al3 { kl ^ k2 ) , 

®ap(ro\ki,k2,u>) = ^ a /3(r \ki,k 2 ,uj) - <f> a p( ■ 

For translationally invariant tensors in homogeneous and 
incompressible turbulence one can write: 



<J a p{k) = P a p(k)a(k) , 
<j> a p{k) = P a p(k)<f>(k) , 

where P a p(k) is the transverse projector, 

k a kp 



(B3) 



P a p{k) = 5 a(: 



k 2 



(B4) 



It is known jl9| that er(fc) (which is the mass operator 
taken at uj — 0) is purely imaginary 



a(k) = —ij(k) 



(B5) 



with 7(fc) real positive. On the other hand 4>{k) is purely 
real. The diagrammatic series expansion of both j(k) 
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and 4>{k) converge order by order, and using scaling rela- 
tions as shown in (|l^) one can find their scaling behavior. 
The order-by-order theory dictates a K41 evaluation of 
these objects which is 



j(k) = c y [ek] 2 ^ , 
4>(k) = c^eAT 3 , 



(B6) 
(B7) 



where c 7 and are dimensionless constants. 

The Dyson- Wyld equations can be written shortly as 



{u)+ivk?)G = P + £*G 
F = G*{& + D)*G, 



(B8) 
(B9) 



where v is the molecular viscosity, P is the transverse 
projector, and D is the correlation function of the ex- 
ternal force which is localized in the energy containing 
interval. The symbol * stands for summation over tensor 
indices and integration over intermediate k. Substituting 
S from Eq. (B3) into the Dyson equation we rewrite: 



(BIO) 



[uj + ivk 2 + ij{k)]G = P + £ * G 



In the bulk of the inertial interval we can neglect vk 2 
with impunity. The zero order solution of this equation 
is obtained by neglecting S: 

G a p -> g a fi(k,u) = P a p{k)g(k,uj) , (Bll) 

9(k,oj)= ] ... ■ (B12) 
uj + i^j(k) 

The zero order solution of F is obtained in three steps: 
first replace <fr by </>, secondly neglect D in the inertial 
interval in comparison with <fi, and lastly substitute g 
instead of G in Eq. (B9). The result is 



F a p -> f a p(k,uj) 



P af ,{k)f(k,w) 

UJ 2 + 7 2 (fc) 



(B13) 
(B14) 



Iterating Eqs. (BE, BE) without the bare forcing and 
viscosity results in a new diagrammatic series, which 
topologically is exactly the same as the old Wyld di- 
agrammatic expansion before line re-summation. The 
difference is twofold. First, instead of bare propagators 
we have K41 propagators g and /, and every 1-particle 
reducible fragment of any diagram will have a counter 
term which subtracts its "diagonal" part. This counter 
term is of no consequence for our procedure here since 
the diagrams involving it are resummed in the 4-point 
vertices (the rungs) together with all the other contribu- 
tions as explained in Appendix |A|. The resulting topolog- 
ical structure of the ladder diagrams is thus unchanged 
in the new formulation. 



APPENDIX C: SELF CONSISTENCY AT THE 
LEVEL OF K41 

Before establishing this self consistency we need to pass 
from correlation functions in k, uj representation to struc- 
ture functions. The theory is done naturally in k, uj rep- 
resentation but the experimental scaling exponents are 
measured in simultaneous structure functions. We first 
transform from ^-representation of p th -order correlation 
function J- p ({kj, ujj}) to simultaneous correlation func- 
tion F p ({kj}) by the integration: 

1 P duj 

F P {{kj})= J H^S(uj 1 + ...uj p )F p ({k j ,uj j }) . 

-oo i=1 

(CI) 

Here {kj,uij} and {kj} are sets of corresponding vari- 
ables with j — 1, . . .p. The transformation from k repre- 
sentation of F p ({kj}) to the p th -order structure function 
is done as follows: define the longitudinal component of 
the velocity as 

S P (r) = {{[-u(^)-u(~)]-^y) . (C2) 

Each of the factors is Fourier transformed according to 



Wi>-«<-5>] -/75W. 



exp ( 



exp ( 



Accordingly, 



S p (r) = (2tt) 



Here 
(2n) s F p ({k j })6(k 1 



■ k p) = ( Y[u( k j) ■ - 

3=1 1 



(C3) 
(C4) 

(C5) 
(C6) 



The functions f p {r, {kj}) are seen from Eq. (C4) to be: 



fp(r,{k 1 }) = Hl^sH^k r r) 

3=1 



In the limit r — > 



f P (r,{k 3 })^l[(krr) 

3=1 



(C7) 



(C8) 



The K41 scaling exponents y p associated with p th -order 
correlation function J- p ({kj, u>j}) oc k~ Vp in (k,u>)- 
representation is 
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y p = Ap- 11/3 



(C9) 



This corresponds to S p (r) oc r p / 3 under the condition of 
convergence of integrals (CI, C5). 

Next consider the 3rd order Green's function, 
G3,3({fcj, Kj}) in which we denoted by kj the set of in- 
coming wave vectors and by Kj the set of outgoing wave 
vectors. The skeleton diagram of G| 3 ({kj , Kj }) which in- 
volves 4-point rungs is shown as diagram (3) in Fig. ||a. 
(The contribution of 6-point rungs to the skeleton is con- 
sidered in Appendix ^ and shown not to change the 
present considerations). This skeleton has two rungs, and 
we consider it in the limit that the incoming kj vectors 
are much larger than the outgoing Kj. In this limit we 
have four Green's functions with large k, contributing 
7jT 4 , and one vertex with all k large, contributing k. The 
two rungs have large k vector in them [k$ in Eq. (p2h], 
giving k 6 . Finally, one of the rungs has large k coming 
and going, and Eq. ( |22] ) requires for it a k 2 ^ 3 . Altogether 
this gives Gs,a({kj , Kj}) oc k X3 with x 3 = 25/3 which 
is equal to 2/3 given by Eq. (C9). This means that the 
skeleton diagrams for G^^({kj 1 Kj}) (with asymptotics of 
the rung defined by the two-point fusion rules) automat- 
ically reproduces the K41 scaling exponent £3 = 1 in the 
three-point fusion. This i s true subject to the condition 
that the integrals (CI, 05) for p = 3 converge. That this 
is so may be shown by a direct calculation. For future 
purposes it is extremely important to note that the prin- 
cipal contribution to the k integral (|C5] ) comes from the 
region where k\ ~ ki ~ £3 ~ 1/r. 

Now let us compare diagram (3) in Fig. ||a and Fig. || 
with the skeleton diagrams for G-i^ and G^. One rec- 
ognizes that in general for G PtP> we will have (p — 1) 
rungs with large incoming k, contributing fc _3 (p _1 ) [orig- 
inating from k e in Eq. (|22|)]. We will have also 2p — 2 
Green's functions with large k contributing fc~( 2 P~ 2 ) 2 / 3 . 
Next we will have p — 2 outgoing legs with large k con- 
tributing /c~(p~ 2 ) 2 / 3 from Eq. (9). Finally we will have 
2p — 2 vertices having incoming and outgoing large k 
vectors, contributing k p ~ 2 . All together we find that 
G„.n({kj, oc k~ Xp with x p = Ap— 11/3 which is equal 



*P,p({kj, K j}) of k 
to y p given by Eq. (C9). Convergence of the k integral 
(C5) for p = 4 may be shown by direct calculations. A 
proof of convergence of the k integrals (C5) for p > 4 is 
a tedious exercise which nevertheless may be done, for 
example, iteratively. It is readily demonstrated that the 
integral converges when all k 2 -vectors are of the same or- 
der of mag nitude (say, k). Then G p . p ~ fc 11 / 3 " 4 ?. After 
{jp — 1) ^-integrations (each of them giving a factor k 2 / 3 ) 
one has £; 3 (p _1 )~p/ 3 which is enough for convergence of 
(p — 1) d 3 k integrals in the UV region kj ~ k ^> 1/r. 
In the IR region kj ~ k -C 1/r the functions f p pr ovid e 
the integral with additional k p factor [according to (JC^)] 
which guarantees the convergence. 

The considerations of the 6-point and higher order 
rungs leave these conclusions invariant. 



APPENDIX D: ANALYSIS OF TWO LOOP 
INTEGRALS CONTRIBUTING TO ( 2 



The integrand in the integral (|6q ) is a function of gi 
and 92 and it depends on k and k as parameters. The 
integration range is the q± — q 2 infinite plane, but in the 
limit k 3> k the main contribution comes form the four 
finite quadrants k < \qi\, \q 2 \ < k. Well inside the quad- 
rants we are allowed to use the asymptotic form in which 
K I <?i 1, 1 92 1 <C k. In this regime the integrand is k, k- 
independent, and the dependence of the integrals on k, k 
appears only via the limits of integration. By changing 
the dummy variables qi and q% we can now project all 
four quadrants into one of them, say q\ and qi positive. 
In this asymptotic regime we can use for the rungs in the 
integrand o f (p6| ) that include either k or k their asymp- 
totic form (|32[). This results in 

K(k,K) = ~S 2 I" ^(11,12) , (Dl) 

Jp 9i Jp 12 

= *(<h,«3)-*(-9i,«a) • (D2) 

In Appendix |e] we show how to analyze this kind of in- 
tegral with the aim of extracting the coefficients of the 
leading and first subleading logarithmic terms, i.e. 



K\ (k, k) = y 1u ' ( k / K ) + b i M*A0 



(D3) 



Using the results there (E4) with * ((71,(72) 

*l(?l,?2), 



*i(oi,t?2) = 



<7 2 |sign((72) 



2(<7f - <?i<?2 



4) 2 



(D4) 



one find immediately a — > ai — 1 as required by the 
anticipated expansion employed in Eqs. (^J)-(^). To 
compute b\ we examine the integral b\(A) numerically, 
see Fig. We see that the requested limit exists and 
that 61 w -0.434. 



X 



-0.41 




-0.42 



-0.43 



FIG. 7. Numerically computed dependence of bi(A) with 
A — exp(a;). 
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APPENDIX E: EXTRACTION OF THE 
SUBLEADING LOGARITHMIC TERM FROM 
THE 2-LOOP INTEGRALS 

The 2-loop integrals have the characteristic structure 



appearing in (Dl) 



Ji 9i Ji 92 



(El) 



where A 3> 1 and ^(gi, 92) is homogeneous function of 
degree zero: *(Aqi,Aq 2 ) = * (91,92)- When * (91,92) = 
1 then 1(A) = In 2 A. In general only the leading term of 
1(A) is proportional to In 2 A and we expect the following 
subleading terms: 



1(A) = (5 2 [-ln 2 A + 61n^ + c+ 4 
2 A 



(E2) 



Our goal is to find the coefficient b in the limit A — > 00. 
Taking the first derivative of (E2) with respect to A and 
multiplying by A we find 



a In A + b 



A 



dq 2 
92 



*(1,92) 



T^(9i,l) 
Ji 9i 



Ji/A x Ji/a y 



(E3) 



where we changed the dummy variables 91 — xA and 
92 = yA. Taking another derivative and multiplying by 
A we find for large A 



Substituting this result in (E3), and representing In A as 
fi/A dx/x we find 



b = lim b(A) , 

A^oo 



dx r 



(E5) 



1 



1 



b(A)= / — *(x,l) + *(M)-tf( 7 ,l)-*(l )7 ) 



If the expansion assumed in Eq. (E2) is valid, this limit 
must exist. 



APPENDIX F: THE NINE 2-LOOP DIAGRAMS 

OF G 3 , 3 

Consider diagram (la) in Fig. M. We are interested in 

(21 (21 

the ratio of T3 la /T|, where T3 / a is obtained by substi- 
tuting the diagram (la) instead of (^3,3 in Eq. ([l9|). In 
the asymptotic regime C 1 the loop integrals over 91 
and 92 contribute mostly in the regime k 3> 91, 92 S> 
In this regime the integrals over k a , w a , fc& , cob cancel in 
the desired ratio. Similarly the Green's functions G(k c i), 
G(K e ) and G(nf) also cancel in the ratio. Accordingly 



(la) 



AAA/ AAAAAAAA/WV- 



AAA/ 



AAA/ AAA/ AAA/ T>WV 



c AAAAAA/- 



1 2 
_AAA UAAA/ 



JAA/V 



(lb) 



* AAA/ 

^ AAA/ 

k, 



AAAAAA^ 
AAA/ AAA/ _ 



^ c AAAAAA/- 



q 2 r 

AAA/ 

-q 2 

AAA/ L 



_AAAAAAAAA- 



AAA 

AAA/ 



(lc) 



AAA/ 

^ AAA 

kc AAAAAA/- 



"AAAAAAa- 



— -^AAA/ TTVW 



AAA 



AAAAAAA/^ ^ • 



AAA/ L/WV LJAAA/ 1 



FIG. 8. The first group of three 2-loop diagrams appearing 
in the loop expansion of Gz^. 



( 2) 

T 3 i a /TJ can be calculated from the amputated diagram 
(E4) (2) in Fig. [l^, in which the explicit dependence on kj and 
u has disappeared. These wavevectors remain only in 



the limits of the integrals over 91 and 92, with k replaced 
by 1/r. In this diagram every black dot contributes a 



factor of 9 



l/3+5„ 



where qj is the wavevector on the right 



of the black dot. This is a remnant of the corresponding 
rung before the amputation. The thin line connecting 
these dots is just a reminder that we have loop integrals 
to perform. 

The point to understand now is that if we use diagrams 
(lb) and (lc) in Fig. | to form T 3 (2) fc and T 3 (2) c , the ratio 
of these to T| can be again calculated from the amputa- 
tion of their own diagrams. This will lead to the identical 
amputated diagram (2) of Fig. 10. In addition, and most 
importantly, the integral that needs to be computed is 
the same as Eq. (|66|). Thus one recaptures Eq. ( |67|) but 
with the combinatorial factor 3 in front of the RHS: 



T. 



(2) 



3,la+lfc+ 



lc( r > K ) 



3J 2 [ ] 
61 In 



In" 
1 

KT 



1 

nr / 
TI(v,k) 



(Fl) 



with b\ of Eq. (pq). 
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* AAA/ 

^ AAA/ 



(2a) 



AAAAAAAAA- 



q 2 



q i r 

UAAA/ AAA/ 



c AAAAAA/- 



AAA/ 



-(q i+q 2 ) 

_AAA/ U 



AAA/ 

_A/W 



AAAAAA/- 



(2b) 



* AAA 

^ AAA/ 

k 



(q i+q 2 ) n K d 

AAAAAA^ 



AAAAAA/ 

-k qi 

JAAAr-^^WV UWV — — 

q 2 



°* a/wwv- 



JAAAAAA^- 



AAA/ 



JA/W — 




c AAAAAAy 



AAA/ 



b AAA 



(3a) 



AAAAAAAAA- 



q 2 



AAA/ 



UAAA/ AAA/ AAAAAA/- 



^ c AAAAAA/- 



-1i 

_AAA/ U 



-q 2 

AAA/ — — 



UAAA/ 



(3b) 



1 AAA/ 



^ b AAA/ 



AAAAAAa- 



qi 



-k c -(q i+q 2 ) 

□AAA/ AAA/ UA/VWW 



c AAAAAA/- 



-q 7 . n 



K 



d 



AAA/ 1 AAA/ 

K 



JAAAAAA/- 



JaaA/ — 



(3c) 



* AAA/ 

^ AAA/ 

k 



AAAAAAa- 



AAAAAA/- 



AAA/ T>WWW- 



^ c AAAAAA/- 



q2 



-(q i+q 2 ) 

AAA/ 



-q 2 



AAA/ 



AAA/ AAA/ 



FIG. 9. The second group of six 2-loop diagrams appearing 
in the loop expansion of G3,3. 



1 ,_, 

■ww — r " 



2 



+ 



91 r-, ^2 

'AAA/ AAA 

-9 1 2 

AA/V- — I 



+... 



+ p(p-D(p-2) 




\ ^AAAAAAr 



FIG. 10. The amputated diagrams that appear up to 
2-loop in the loop expansion of G PtP , with the appropriate 
combinatorial factors. 

The second group of six diagrams (2a - 3c) shown in 
Fig. [)] yields to a similar analysis, but the amputated dia- 
gram is shown as diagram (3) in Fig. [HI All six diagrams 
result in the very same amputation, up to permutations 
of the three struts. Analyzing the amputated diagram (3) 
one brings it to the canonical form ( D 1 ) with \f (qi , q — 2 ) 
given by (|70[). Accordingly we write 



lS.,...,>,«) = 6J 2 



(") 

V ft?'/ 



V 

2 \nr * 



6 9 In — 



(F2) 



with 6 2 of Q7ip. 

The analysis of the 2-loop diagrams that involve 4- 
point rungs in the context of G p . p follows exactly the 
same lines, with the amputated diagrams being those of 
Fig. [l^. The only thing to mind is the combinatorics, 
which are presented explicitly in Fig. 10, leading to the 
numbers in Eq. ([7q). 



APPENDIX G: RESUMMED EQUATIONS FOR 
THE 4-POINT AND 6-POINT RUNGS 



In this Appendix we sketch a theory for the 4-point 
and 6-point rungs. Our main aim here is to explain why 
the 6-point rung is quadratic in the smallness, but we use 
the opportunity to indicate how a future theory of these 
objects may be formulated. 

Consider the beginning of the series expansion of the 
4-point rung which is shown in Fig. ^b. Diagram (2) con- 
tains a cross of correlators each attached to two 3-point 
vertices. This is exactly diagram (1), and therefore the 
equation can lend itself to resummation resulting in the 
equation shown in Fig. pd]a. We note that this is not the 
full equation for the 4-point rung even in 1-loop order 
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FIG. 11. A sketch of the equations for the 4-point rung 
(panel a) and 6-point rung (panel b). The shading of rungs 
is used to identify different pieces of the same object. In the 
same manner we shade the loops to identify one or two loop 
integrals. 



since we did not take into account the ladders with a cor- 
relator and Green's function in a cross section. Taking 
into account all the needed contributions is not difficult, 
but is not the main point of this Appendix, and we pro- 
ceed for simplicity without the additional terms. 

In the asymptotic regime the bare contribution dia- 
gram (1) in Fig. pr| a is negligible. With this contribution 
discarded, the remaining equation is homogeneous, call- 
ing for finding a zero mode of the equation. Since we have 
already demonstrated that the 4-point rung is small, of 
the order of 62 we can conclude that the loop integral 
which we denote as £\ must be large, or the order of 
1/(52- (the homogeneous equation can be only solved if 
<5 2 Ri b\l\). In fact, in the future it would be extremely 
worthwhile to solve the full equation in the 1-loop or- 
der and demonstrate that this is the case, and thus to 
lend further weight to the theory presented in this pa- 
per. Of course solving such an equation will also supply 
us with a functional form of the 4-point rung, and with it 
a substantial part of the value of the parameter 62 which 
appears in the final result for the scaling exponents. 

In Fig. pT) p we present the resummed form of the equa- 
tion of the 6-point rung, to the same level of qualitative 
discussion. Again we discard in the asymptotic limit the 
bare contribution of diagram (1), but we cannot neglect 
diagram (2) since it has the same asymptotic behavior as 



the resummed 6-point rung. Diagram (2) is of the order 
of b\ . Diagram (3) is of the order of ec>2^2 where £2 is 
the loop integral. This integral is very similar to t\, and 
we therefore estimate £2 rs i\ ~ 1/5 2 t and thus diagram 
(3) is of the order of the LHS. Diagram (4) is of the or- 
der of e 2 i?3^4 where £3 and £4 each refers to one of the 
loop integrals. With the same level of approximation we 
estimate it thus to be of 0(e 2 /<5 2 ). Denoting x = e/S 2 
we thus represent the order of magnitude relations that 
result from panel b by the equation 

x = 1 + ax + bx 2 , (Gl) 

where a and b are dimensionless constants of 0(1). It is 
obvious that only x ~ 1 is a consistent solution of this 
equation, and we thus conclude that the 6-point rung is 
quadratic in the smallness 62- 

We therefore understand that the 6-point rung appears 
in our considerations only at the level of the 0(S 2 ) order. 
In this order it appears in addition to the 2-loops integrals 
which are formed by two 4-point rungs, as discussed in 
detail in the text of the paper. But since the 6-point rung 
connects three struts, exactly like the structure made 
of two 4-point rungs, the combinatorical factors appear- 
ing in the pth order scaling exponents are identically the 
same. Accordingly we understand that the effect of the 
6-point rung is only in renormalizing the value of the 
parameter 62 which anyway is model dependent. 

Similar consideration apply to the 8-point rung which 
begins to affect the theory only in 0(<5|). It will renor- 
malize the value of the parameter 63 in Eq. ( |78| ) . Higher 
order rungs are even less relevant for the calculation at 
hand. 
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